/3-decay half-lives and /3-delayed neutron emission probabilities of nuclei in the region 

A < 110, relevant for the r-process 
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Measurements of the /?-decay properties of A < 110 r-process nuclei have been completed at the 
National Superconducting Cyclotron Laboratory, at Michigan State University. /3-decay half-lives 
for i^^Y, 106.10721. and ^"Mo, along with /3-delayed neutron emission probabilities of ^'^^Y , lo^.iiOjyjQ 
and upper limits for ^"^Y, lO^-iOT^j. ^j^^j i^^'^Mo have been measured for the first time. Studies on 
the basis of the quasi-random phase approximation are used to analyze the ground-state deformation 
of these nuclei. 

PACS numbers: 23.40.-s; 21.10.Tg; 27.60.+j; 26.30.-k 
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I. INTRODUCTION 



The rapid neutron-capture process (r-process) [l|, Q 
remains as one of the most exciting and challenging 
questions in nuclear astrophysics. In particular, the 
theoretical quest to explain the production of r-process 
isotopes and the astrophysical scenario where this pro- 
cess occurs have not yet been satisfactorily solved (for 
a general review see, for instance, Ref. 0, 0, Q). R- 
process abundance distributions are typically deduced 
by subtracting the calculated s- and p-process contri- 
butions from the observed solar system abundances. 
Furthermore, elemental abundances originated in the 
early Galaxy can be directly observed in metal-poor, 
r-process-enriched stars (MPRES) (i.e., [Fe/ H]<—2, 
[Ba/Eu]<-0.7, [£'u/i^e]>-|-l). These combined obser- 
vations reveal disparate behavior for light and heavy nu- 
clei: MPRES abundance-patterns are nearly consistent 
from star to star and with the relative solar system r- 
process abundances for the heavier neutron-capture ele- 
ments A>130 (Ba and above), suggesting a rather ro- 
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bust main r-process operating over the history of the 
Galaxy. Such a consistent picture is not seen for light 
neutron-capture elements in the range 39<Z<50, as 
the solar system Eu-normalized elemental abundances 
in MPRES show a scattered pattern [1, 0, H, 111 ■ Anti- 
correlation trends between elemental abundances and Eu 
richness at different metalicities have been suggested to 
provide a hint for an additional source of isotopes be- 
low A-UO [13, UU, El, d, [3. Measured abundances 
of i°^Pd and 1291, for A<130, and ^^^Hf, for A>130, 
trapped in meteorites in the early solar system forma- 
tion [11,13, further reinforce the idea of different ori- 
gins for isotopes lighter and heavier than A=130 (see 
e.g. Ref. [13). 

Reliable nuclear physics properties for the extremely 
neutron-rich nuclei along the r-process path are needed 
to interpret the observational data in the framework of 
proposed astrophysical models. The r-process abundance 
region around A~110, prior to the A=130 peak, is one 
region of intense interest, where the astrophysical models 
underestimate the abundances by an order of magnitude 
or more. It has been shown that the underproduction of 
abundances can be largely corrected under the assump- 
tion of a reduction of the N—82 shell gap far from sta- 
bility 19, 20, 21], [22|. Whereas such quenching effect 
was suggested in elements near Sn — via analysis of the 
P decay of ^^^Cd into -'^^^In [3 — its extension towards 
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lighter isotones 47^Ag82, 46*Pd82 and below would have 
crucial consequences in the search for the r-process site. 
In this sense, self-consistent mean-field model calcula- 
tions predict that the N=82 shell quenching might be 
associated with the emergence of a harmonic-oscillator- 
like doubly semi-magic nucleus \}fZT^o, arising from the 
weakening of the energy potential surface due to neutron 
skins 'M, in, HE '27]. Thus, it is imperative to char- 
acterize the evolution of nuclear shapes in the region of 
"OZr. 

The goal of the experiment reported here was to use 
measured /3-decay properties of nuclei in the neighbor- 
hood of ^^°Zr to investigate its possible spherical char- 
acter arising from new semi-magic numbers [l^, [1^ , or 
even a more exotic tetrahedral-symmetry type predicted 
by some authors jll, |2^. In particular, the measured 
half-lives (7i/2) and /3-delayed neutron-emission proba- 
bilities (Pji) can be used as first probes of the structure 
of the /3-decay daughter nuclei in this mass region, where 
more detailed spectroscopy is prohibitive owing to the 
low production rates at present radioactive beam facil- 
ities. A similar app roach has already been used in the 
past [M,!!! SIM [missile. Besides the nuclear- 
structure interest, our measurements also serve as impor- 
tant direct inputs in r-process model calculations. The 
Ti/2 values of r-process waiting-point nuclei determine 
the pre freeze-out isotopic abundances and the speed of 
the process towards heavier elements. The Pn values of 
r-process isobaric nuclei define the decay path towards 
stability during freeze-out, and provide a source of late- 
time neutrons. 

In the present paper, the measurements of T1/2 and 
P„ of 100-105Y, i03-i072r, i06-i09j^b, 108-iiiMo and 
i09-ii3rpp are reported. The work is contextualized amid 
a series of /3-decay r-process experimental campaigns, 
carried out at the National Superconducting Cyclotron 
Laboratory (NSCL) at Michigan State University (MSU). 
Details of the experiment setup and measurement tech- 
niques are provided in Sec. |TT1 followed by a description 
of the data analysis in Sec. IIIII The results are further 
discussed in Sec. |IV] on the basis of the quasi-random 
phase approximation (QRPA) [H, [H, H^, |4l| using nu- 
clear shapes derived from the finite-range droplet model 
(FRDM) [42] and the latest version of the finite-range 
liquid-drop model (FRLDM) [Ulil]. The paper is closed 
with the presentation of the main conclusions and future 
plans motivated by the current measurements. 



II. EXPERIMENT 
A. Production and separation of nuclei 

Neutron-rich Y, Zr, Nb, Mo and Tc isotopes were pro- 
duced by fragmentation of a 120 MeV/u ^^^Xe beam in 
a 1242 mg/cm^ Be target. The primary beam was pro- 
duced at the NSCL Coupled Cyclotrons [45| at an aver- 
age intensity of 1.5 pnA (^1 x 10^° s~^). Forward-emitted 



fragmentation reaction products were separated in-flight 
with the A1900 fragment separator [i^ — operated in its 
achromatic mode — using the Bp-AE-Bp technique [47j . 
Two plastic scintillators located at the intermediate (dis- 
persive) focal plane and at the experimental area were 
used to measure the time-of-flight (ToF), related to the 
velocity of the transmitted nuclei. The first of these 
scintillators provided also the transversal positions Xd of 
the transmitted nuclei at the dispersive plane. A Kap- 
ton wedge was mounted behind this detector to keep 
the achromatism of the A1900. Energy-losses experi- 
enced by nuclei passing through this degrader system 
of 22.51 mg/cm2 (Kapton) and 22.22 mg/cm^ (BC400) 
thickness, provided a further filter to select a narrower 
group of elements. A total of 29 neutron-rich iso- 
topes (100-105Y, i02-i072r, i04-i09Nb, 106-iiiMo and 
i09-ii3rp^-j defined the cocktail beam that was transmit- 
ted to the implantation station. 



The trajectory followed by the nuclei at the A1900 dis- 
persive focal plane depended on their magnetic rigidities 
Bp, which, in turn, were related to the corresponding ve- 
locities and mass-over-charge ratios. An event- by- event 
separation of the transmitted nuclei according to their 
mass A and proton Z numbers was achieved by com- 
bining the measured ToF and Xd with the energy-loss 
AE in a silicon PIN detector located in the implantation 
station. The latter quantity had to be corrected from its 
velocity dependence (described by the Bethc-Bloch equa- 
tion |48| ) ; the identification of nuclei transmitted through 
the A1900 is shown in Fig. [T] In this figure, the vari- 
able ToF* corresponds to the ToF corrected from its Xd 
dependence, and the AE signal was measured with the 
most upstream PIN detector of the implantation station. 



The maximum 5% Bp acceptance of the A1900 in- 
cluded the nuclei of interest, along with three primary- 
beam charge-states ^^^Xe+^^ (Pp=3.8251 Tm), ^^^Xe+^° 
(Sp=3.9016 Tm) and ^^^Xe+^^ (Bp=3.9812 Tm), with 
particle rates of 8.7x10^ s-\ 3.5x10^ s~i and 80 s-\ 
respectively. As these contaminants could reach the in- 
termediate image plane, resulting in the damage of the 
plastic scintillator, it was necessary to stop them with the 
standard A1900 slits and a 17.4 mm-wide tungsten finger 
located in the first image plane. This slit configuration 
blocked the most intense contaminants with a minor re- 
duction of the Bp acceptance of the fragments of interest. 
Thus, the ^'^^Xe"'"^° charge-state was blocked by the fin- 
ger at the central position of the transversal plane, while 
^■^^Xe+^^ stopped in one of the slits, closed at 39.88 mm 
from the optical axis (in the low- Bp side). The high- 
Bp slit was fully opened so that the fragment of inter- 
est (along with the low-intensity contaminant ^^^Xe"*"*^) 
were transmitted through the second half of the A1900. 
The resulting overall Bp acceptance was about 4%. 
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B. Implantation station 

Nuclei transmitted through the A1900 and beam trans- 
port system were implanted in the NSCL beta counting 
system (BCS) [i^ for subsequent analysis. The BCS con- 
sisted of a stack of four silicon PIN detectors (PINl-4) 
of total thickness 2569 /xm, used to measure the energy- 
loss of the exotic species, followed by a 979 /im-thick 
40 X 40-pixel double-sided silicon strip detector (DSSD) 
wherein implanted nuclei were measured along with their 
subsequent position-correlated (3 decays. Located 9 mm 
downstream of the DSSD, a 988 /xm-thick 16-strip single 
sided Si detector (SSSD) and a 10 mm-thick Ge crystal, 
separated 2 mm from each other, were used to veto par- 
ticles whose energy-loss in the DSSD was similar to that 
left by P decays. The signals from each DSSD and SSSD 
strip were processed by preamplifiers with high- and low- 
gain outputs, each defined over a scale-range equivalent 
to 100 MeV and 3 GeV, respectively. Energy thresh- 
olds for the low-gain signals were set to values around 
300 MeV. High-gain signals from every DSSD and SSSD 
strips were energy-matched in the beginning of the exper- 
iment using a ^^*Th a-source as reference, whereas inde- 
pendent threshold were adjusted using a ^°Sr /?-source. 
Two of the DSSD strips (the front-side 31*^' and the back- 
side 12"^) were damaged and had consequently to be dis- 
abled during the experiment, raising their high/low-gain 
thresholds to maximum values. Energy thresholds for the 
PIN detector were set to 2000 MeV for the first detec- 
tor and around 300 MeV for the rest of PIN detectors. 
A dedicated un-wedged setting of the A1900, transmit- 
ting a large amount of light fragments, was used in the 
beginning of the experiment to energy-calibrate the PIN 
detectors and low-gain signals from DSSD. The measured 
energy losses were compared with calculations performed 
with the LISE program [50] , using the Ziegler energy-loss 
formulation [5l| . 

The BCS was surrounded by the neutron emission ra- 
tio observer (NERO) detector [H, |H, [13], which was 
used to measure /3-delayed neutrons in coincidence with 
the /3-decay precursor. This detector consisted of sixty 
(16 ^He and 44 B3F) proportional gas-counter tubes, em- 
bedded in a 60 X 60 X 80 cm^ polyethylene moderator ma- 
trix. The detectors were set parallel to the beam di- 
rection and arranged in three concentric rings around 
a 22.4 cm-diameter vacuum beam line that accommo- 
dated the BCS. /3-delayed neutrons were thermalized in 
the polyethylene moderator to maximize the neuron cap- 
ture cross section in the "^He and B3F gas counters. Gains 
and thresholds of the neutron counters were adjusted us- 
ing a ^^^Cf post-fission neutron source. 

The master-event trigger was provided by the PINl de- 
tector or by a coincidence between the DSSD front and 
back high-gain outputs. The master trigger opened a 
200 /xs-time window [13, [13l , during which signals from 
the NERO counters were recorded by an 80 ms-range 
multi-hit TDC. The 200 /xs interval was chosen on the 
basis of an average moderation time of about 150 /xs. 



measured for neutrons emitted from a ^^^Cf source. En- 
ergy signals of moderated neutrons were also collected in 
an ADC. The closure of the time window was followed 
by the readout of the BCS and NERO electronics. 



C. Identification of nuclei 

The particle identification (PID) was performed with 
a dedicated setup installed upstream of the BCS/NERO 
end station. Several nuclei having /xs isomers with known 
7-decay transitions were selected in a specific A1900 set- 
ting, and implanted in a 4 mm-thick aluminum foil, sur- 
rounded by three 7-ray detectors from the MSU seg- 
mented germanium array (SeGA) [1^. These detectors 
were energy and efficiency calibrated using sources of 
^^Co, ^°Co and ^^^Eu. The 7-peak efficiency was about 
6% at 1 MeV. A 50x50 mm^, 503 /xm-thick sihcon PIN 
detector (PINO) upstream of the aluminium foil was used 
to measure Ai5, which, combined with ToF and Xd, al- 
lowed for an identification of the various nuclei transmit- 
ted in the setting according to their Z and A numbers. 
The measurement of these three signals, gated in specific 
known isomeric 7-decay lines, provided a filter to identify 
the corresponding /is isomers. 

Four settings of the A1900 were used to identify the nu- 
clei of interest in a stepwise fashion based on the observa- 
tion of the ^s isomers. In each of these settings, the first 
half of the A1900 was tuned to 5/9=3.9016 Tm— defined 
by the rigidities of the nuclei of interest — while the sec- 
ond half was tuned to Bp=3.7881 Tm, Bp=3.7929 Tm, 
Bp=3.7976 Tm and Bp=3.8024 Tm. The four val- 
ues were chosen to have fragments transmitted in two 
consecutive settings, and far enough to cover the en- 
tire range of rigidities of interest. In the first setting 
i3p=3.7881 Tm, several 7 lines were seen for the /is- 
isomeric states of ^^ipd (135 keV), ^^SAg (714 keV), 

Ag (156 keVV^i^^Ag (684 keV) and ^^^Cd (720 keV 
and 743 keV) [5^, [13|. From these references, it was 
possible to identify the more exotic nuclei present in the 
subsequent settings. The particle identification was fur- 
ther confirmed by detecting the 135 keV 7 line from the 
isomer ^^^Pd transmitted in the settings i3p=3.7929 Tm 
and Bp=3.7976 Tm. The fourth setting Bp=3.8024 Tm 
was chosen to optimize the transmission of the nuclei of 
interest to the final end station. 

After the PID was completed, the Al catcher and PINO 
detector were retracted. The cocktail beam was then 
transmitted to the final BCS/NERO end station, where 
the Ai? necessary for the PID spectrum was provided by 
the PINl detector. 

The PID spectrum shown in Fig. [1] includes the fully- 
stripped ions of the nuclei of interest, overlapped with a 
small fraction of charge-states contaminants from lighter 
isotopes. Given the high Bp selected in the A1900 set- 
ting, only the hydrogen- like ions with mass numbers A— 2 
and j4 — 3 reached the experimental area. These contami- 
nants were disentangled from the fully-stripped nuclei by 
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ToF (arb. units) 

FIG. 1: (Color online) Particle-identification (PID) spectra 
collected during 112 hours of beam time, showing the different 
nuclei characterizing the cocktail beam transmitted to the end 
station. The solid line marks the r-process region of interest 
(left side) . The r-process waiting-point nuclei analyzed in the 
present experiment are indicated by the ellipses. 



measuring the total kinetic energy (TKE) of the trans- 
mitted ions compiled from signals of the PIN detectors 
and DSSD. The TKE spectra of different Zr isotopes is 
shown in Fig. O the upper row of panels corresponds to 
the nuclei that were detected with the first PIN detec- 
tor of the BCS. The double-peak structure in the TKE 
spectra arises from the fully-stripped species {In^-TKE 
peak) and the corresponding hydrogen- like contaminants 
(low-TKE peak). The central and lower row of panels 
show the same spectra with the additional requirement 
of being implanted in either the most downstream PIN 
detector in the stack (i.e., PIN4) (central row) or in the 
DSSD (lower row). Nearly only the fully-stripped nuclei 
reached the latter, whereas the hydrogen-like components 
were mainly implanted in the last PIN. A small fraction 
of the fully-stripped components did not reach the DSSD 
due to a slightly overestimated Si thickness. Further- 
more, no low-gain signals from the SSSD were observed, 
demonstrating that no nuclei reached this detector. 



III. DATA ANALYSIS 

A. /3-decay half-lives 

Specific conditions in the different detectors of the BCS 
were required to distinguish implantations, decays and 
light-particle events: a signal registered in each of the 
four PIN detectors, in coincidence with the low-gain out- 



put from at least one strip on each side of the DSSD, and 
in anti-coincidence with the SSSD, was identified as an 
implantation event. Decay events were defined as high- 
gain output signals from at least one strip on each side 
of the DSSD, in anti-coincidence with signals from PINl. 
Software thresholds were set separately for each DSSD 
strip in order to cut off noise. Decay-like events accompa- 
nied by a preamplifier overflow signal from the Ge crystal 
downstream of the DSSD were identified as light parti- 
cles and consequently rejected. According to LISE-based 
calculations, these light particles were mainly tritium nu- 
clei and, to a lesser degree, ^Li nuclei, with energy-loss 
signals in the DSSD high-gain output comparable to the 
decays of interest. 

For each implantation event, the strip location on each 
side of the DSSD — defining the implanted pixel — was 
determined from the average of the strips weighted by 
their respective energy signal amplitude. The resulting 
average pixel was recorded along with the implantation 
time taken from a continuously counting 50 MHz clock. 
The last beam-line quadrupoles in front of the BCS were 
adjusted to illuminate a wide area of the DSSD cross- 
section; the resulting distribution of implantation events 
is shown in Fig. [31 Subsequent decay events occurring 
in the same or neighboring pixels (defining a cluster of 
nine pixels) within a give correlation-time window (tc) 
were associated with the previous implantation, and their 
times and pixels recorded. The value of tc was chosen 
to be around 10 times the expected Ti/2- Whenever 
a decay-like event was correlated with more than one 
implantation, all the events within the sequence (i.e., 
decay and implantations) were rejected. Such scenario 
would be possible if different implantations occur in the 
same cluster within the correlation-time window. Given 
the low maximum implantation rate per pixel of about 
1.8 xlO^'^ s^^ and the T1/2 values (typically below 1 s), 
the probability of multiple implantation events correlated 
with a decay was negligible. 



1. 13 -decay background 

The implantation-decay correlation criterion did not 
prevent the occurrence of spurious correlations arising 
from sources other than the actual decays of interest. 
Possible background sources included: light particles 
that did not lead to overflows in the Ge detector (either 
because they missed the detector or because they de- 
posited only a fraction of their energy) ; real decays from 
longer-lived implanted nuclei and from nuclei implanted 
in neighboring pixels; and electronic noise signals above 
the thresholds. A detailed study of this decay- like back- 
ground was necessary to extract T1/2 values. 

Background rates were determined separately for each 
DSSD cluster and 1-hour data-collection run, counting 
the number of decay-like events that were not correlated 
with any implantation. During this background measure- 
ment, each time an implantation was detected in a given 
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FIG. 2: Total kinetic energy (TKE) of Zr isotopes identified in the PID spectrum. The different spectra rows correspond to: 
aU detected events in PINf (upper row); nuclei that were implanted in PIN4 (central row); and nuclei that were implanted in 
DSSD (lower row). The fully-stripped ions are labeled on top of the corresponding high-T KE peaks. 




time (i.e., those that were not correlated with the previ- 
ous implantation) were recorded as background events. 
Background rates were calculated for each 1-hour run as 
the ratio of the number of uncorrelated decay events in 
each cluster to the unblocked time in that cluster. The 
resulting rates were position-dependent and nearly con- 
stant over different runs. A critical factor in this analy- 
sis was the length of the post-implantation blocking time 
which had to be chosen so as to minimize the proba- 
bility of recording real correlated decays. A blocking- 
time window of 40 s was found to fulfill this requirement. 
Finally, the DSSD cluster-averaged /3-decay background 
was about 0.01 s~^, nearly constant throughout the ex- 
periment. 



FIG. 3: (Color online) Implantation-event distribution mea- 
sured over 112 hours of beam on target. The damaged 12**^ 
back strip and 31°' front strip in the DSSD can be observed 
as respective gaps on the top-plane projection. 



pixel, the corresponding 9-pixel cluster was blocked to 
any subsequent decay event during a time interval chosen 
to be longer than tc- Decay-like events detected in that 
cluster after the closure of the post-implantation blocking 



The total number of /3-decay background events Bp 
for each isotope was calculated from the measured back- 
ground rates in the runs and pixels where the isotopes 
were implanted, multiplied by tc- This number was di- 
vided by the total unblocked time, calculated as the prod- 
uct of the number of implantations and tc, giving the 
background rate for that nucleus. The statistical error 
for the background rate of each nucleus, derived from 
the number of background events recorded in each DSSD 
cluster and run, was about 5%. 
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2. 13-decay half-lives from decay-curve fits 

The time differences between impiantation and corre- 
iated decay events were accumulated for eacfi nucleus 
in separated histograms and fitted by least-squares to 
a multi-parameter function derived from the Batemann 
equations ^5^. Due to the limited detection efficiency 
of the DSSD, some of the recorded /3-decay events may 
come from descendant nuclei following a missed decay 
of the nuclei of interest. Given the low probability of 
missing consecutive decays, the correlation times used in 
the analysis, and the values of the half-lives of the de- 
scendant nuclei, up to three generations were included 
in the fit functions, along with contributions from back- 
ground events. The paths that define the possible decay 
sequences following the decay of a mother nucleus are 
schematically illustrated in Fig. 21 
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FIG. 4: Schematic representation of the possible decay-paths 
followed by a decaying mother nucleus (1), including first (2) 
and second (3) descendant-nuclei generations (see text for de- 
tails). 



FIG. 5: Decay curve data (full circles) for ^''''■^''^Zr64,65 iso- 
topes. Included are fit functions for mother (solid thin line), 
daughter (dashed line), granddaughter (dot-dashed line), 
background (dotted line). The sum of these functions is rep- 
resented by the solid thick line. 



The fit equation included a total of eleven parameters, 
eight of them fixed to constant values, namely: the decay 
constants of the daughter (A2) and granddaughter (A3); 
the neutron-emission probability of the mother (-Pin), 
daughter (p2n), and neutron-emitted daughter (P2nn); 
and the decay constants of the single neutron-emitted 
daughter (A2n) and granddaughter (Asn), and the double 
neutron-emitted granddaughter (Asnn)- These fit con- 
stants were taken from the literature or, in the case of 
some unknown Pin, calculated using the FRDM+QRPA 
model [3^, HO, Si], . The remaining three parameters 
were treated as free variables to be determined from the 
fit algorithm; two of them were the decay constant of the 
mother nucleus (Ai) and the initial number of mother 
decaying nuclei (A^o)- The third parameter, namely the 
background constant, was treated as a constrained "free" 
variable, defined within ±10% of the calculated value 
(see Sec. IIII A II) . Finally, after determining the Pn as 
described in Sec. IIIIBi the decay curves were re-fitted 
replacing Pin by the newly measured values. The decay 
curves of some selected Zr isotopes are presented in Fig.[S] 
with the different contributions to the total fit curve. 



The least-squares method used to fit the decay-curve 
histograms requires that the number of events per bin size 
Ai is described by Gaussian statistics. More formally, the 
individual probability for one event to be recorded in a 
given bin must be <ICl, and the total number of events per 
bin A(Ai) must be large, typically ~20. Taking the time 
scale of the histograms as tc, the latter condition can be 
expressed as N{At) = {At/tc)N - 20, where A^ is the 
total number of decay events in the histogram. Thus, 
since At/tc -Cl, A must be ^200 for the least-squares 
method to be valid. Table |T] shows the half- lives of those 
nuclei that fulfilled the Gaussian-statistics requirement, 
along with their corresponding A. For cases with lower 
statistics, an alternative analysis based on the Maximum 
Likelihood method was used (see Sec. IIII A^ . 

Different sources of systematic error were included in 
the decay-curve analysis: uncertainties in the input pa- 
rameters (half-lives and neutron-emission probabilities of 
the descendant nuclei) were accounted for by comparing 
the fit results obtained with these input values scanned 
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FIG. 6: /3-decay half- lives obtained from the Maximum Likelihood method for Y, Zr, Nb, Mo and Tc isotopes (filled circles), 
compared with results from previous experiments [3^ . |33| . [sol . IgoI . [61| . [g^ . [g^ . |64, 65, 6^, [gtI . [G^ . [gqI [70| (open circles). For the 
sake of clarity, the latter were shifted to the right by 0.1 units. The data are compared with two versions of the QRPA model 
of MoUer et ai: the version described in Ref. [4l[ (solid line), and the interim version QRPA06 described in Sec. II VI (dashed 
line). (See text for details). 



over their respective error intervals. The resulting errors 
depended on the half-lives of the mother and, to a lesser 
degree, descendant nuclei. Uncertainties were typically 
below 5%. In addition, comparisons of fit half-lives using 
background rates varied over their corresponding uncer- 
tainty showed differences below 1%. Absolute systematic 
and statistical errors are shown in Table [D 



3. /3-decay half-lives from Maximum Likelihood Method 

The Maximum Likelihood analysis (MLH) is well 
suited for determinin g de cay half-lives in cases with low 
statistics [H llllzl H, S [ll]. The method, de- 
scribed in appendix [XI defines decay sequences for up to 
three generations following an implantation. The proba- 
bility of observing a given decay sequence was calculated 
by summing up the probabilities for all possible scenar- 
ios leading to the detection of the decay-event members. 
The scenarios were evaluated by considering the occur- 
rence of up to three (3 decay events, including /3-delayed 
neutron branching, the contributions from background 
events, and the "missing" decays due to the limited detec- 
tion efficiency. A joint probability density, the likelihood 
function C, was calculated by multiplying the probabili- 



ties for all the measured decay sequences. The resulting 
C was a function of the measured decay times U of the 
different members of the decay-sequence, their decay con- 
stants and neutron-emission probabilities, the correlation 
time tc, the background rate of the corresponding DSSD 
cluster and run where the decay sequence was detected, 
and the /3-decay detection efficiency e^j. The half- lives 
of the nuclei of interest were determined from the max- 
imization of C, using the decay constant of the mother 
nucleus Ai as free parameter. 

All of the descendant decay parameters necessary to 
define C were taken from previous measurements or — in 
the case few values — calculated from theory. Similar 
to the least-squares fit method described in Sec. IIII A21 
the half-lives were re-calculated with the newly measured 
Pn values, once known. The /3-decay detection efficiency 
6/3 was determined as the ratio of the number of detected 
/3-decays Nf} attributed to a given nucleus to the number 
of implantations of that nucleus. The former was given 
by = No/{XiAt), where A^o and Ai were obtained 
from the decay-curves fits, for the cases where the least- 
squares method was valid. No systematic trend for 
was observed within a given isotopic chain, so a weighted 
average efficiency per DSSD cluster of (31±4)% was used. 
Finally, the background rate was determined for each 
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TABLE I; Total number of implantations, number of events per histogram A*" and experimental /3-decay half-lives obtained 
from least-squares fit (Least squares) and the Maximum Likelihood method (MLH) with systematic and statistical errors. The 
results are compared with available data from previous experiments (Literature) , and with the versions QRPA03 and QRPA06 
of MoUer's QRPA model. (See text for more details). 
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DSSD cluster and run, as described in Sec. IIII A II 

The sources of systematic error included contributions 
from uncertainties in the experimental descendant-nuclei 
Ti/2 and Pn, background and ep. The systematic error 
of T i/2 was calculated for each nucleus as described in 
Sec. IIII A2[ yielding typical values below 10%. The sta- 
tistical error was directly calculated from the MLH anal- 
ysis using the prescription described by W. Briichle [77j . 
Since the C distributions were typically asymmetric, the 
shortest possible interval containing the maximum of the 
/^-distribution and 68% (i.e., 1-a) of the total integrated 
density probability was used [13, ■ The calculated sys- 
tematic and statistical uncertainties are listed in Table [J 
The total error shown in Fig. [5] was obtained summing up 
the contributions from systematic and statistical uncer- 
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tainties according to the method described in Ref. fl^ - 
The Ti/2 obtained from the MLH are in agreement with 
the decay-curve fits for the cases where the least-squares 
fit method was valid. 



B. /3-deIayed neutron emission probabilities 

The f] decay of a neutron-rich nucleus can populate 
levels in the daughter nucleus above the neutron sepa- 
ration energy Sn, thus opening the /3-delayed neutron- 
emission channel. The probability of observing a neu- 
tron associated with the f3 decay of a nucleus is given 
by the neutron-emission probability or Pn value (called 
Pin in Sec. IIII A~2|) . /3-delayed neutrons were detected in 
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coincidence with /? decays using the NERO detector in 
conjunction with the BCS. values were determined for 
each nucleus according to: 



/3/3n 



enNp 



(1) 



where Np^ is the number of detected neutrons in coinci- 
dence with P decays correlated with previous implanta- 
tions; i?n is the number of background /3-neutron coin- 
cidences; En is the neutron detection efficiency; and Np 
is the number of /3-decaying mother nuclei. Npp^ is the 
number of detected /3-delayed neutrons from descendant 
nuclei; it should be consequently subtracted from the to- 
tal number of detected neutrons in order to determine the 
actual number of neutrons associated with the nucleus of 
interest. For the nuclear species discussed in this pa- 
per, /3- neutron coincidences associated with descendant 
nuclei other than the /3-decay daughter were negligible. 
Using the Batemann equations (ssj . it is possible to write 
explicitly the value of Npp-c^: 



Nppn = (1 - Pn)C, 



where C is a constant given by: 



C 



A2-A1 



Ai 
A2 



(2) 



(3) 



In this equation, P^n is the neutron-emission probability 
of the daughter nucleus (called P2n in Sec. IIII A2l) . and 
A2 and Ai are the decay constants of the daughter and 
mother nuclei, the latter being extracted from the anal- 
ysis discussed in the previous sections. Inserting Eq. [5] 
and Eq. [3] into Eq. [l] and rearranging terms: 



Pn = 



Nsn — — C 



N. 



c 



(4) 



The value of iVg for a given nucleus was calculated as 
the product of the total number of implantations by the 
average e/3. The number of neutrons detected by NERO 
in coincidence with /3 decays were recorded in a multi-hit 
TDC. Pn values were first determined for the less ex- 
otic nuclei, taking Ai from Table |T1 and A2 and P^n from 
Ref. [6lj (as their corresponding daughters were not in- 
cluded in the present experiment). The newly calculated 
values of Pn were then included in Eq. [?]as Pnn, to cal- 
culate Pn for the next exotic nuclei. 



1. Neutron detection efficiency 

The design of the NERO detector was optimized to 
achieve a large and energy-independent efficiency, at least 
in the typical range of energies of the measured /9-delayed 
neutrons. The efficiency response of NERO was deter- 
mined at the Nuclear Structure Laboratory, at the Uni- 
versity of Notre Dame, by detecting neutrons produced at 




o 
c 
9 
o 
3= 
LU 



10 

E„ (MeV) 

FIG. 7: MCNP-calculated efficiencies as a function of neutron 
energy En for the innermost ring (solid line) , intermediate ring 
(dashed line) and external ring (dotted line), compared to the 
measured values from the reactions ^^]i{a,n) (dots), ^^C{a,n) 
(triangles) and ^^'V{p,n) (empty circles), and from a neutron 
post-fission ^'^^Cf source (square). See Ref. [s^. [Hsj] for more 
details. 



different energies from resonant and non-resonant re- 
actions, and from a ^^^Cf source, as described in Ref. [53 |. 
In that analysis, eight different values of En ranging from 
about 0.2 MeV to 5 MeV were covered. The experimental 
results were extrapolated to a wider energy range, using 
the MCNP code The NERO efficiency is nearly 

constant for En below 0.5 MeV, and gradually decreases 
beyond this value, as discussed in Refs. [H, [5J|. Further 
analysis of the detector rings showed that, for energies 
below 1 MeV, where NERO is most efficient, the total 
efficiency was mainly governed by the innermost detec- 
tor ring followed by the intermediate and external rings. 
This result suggested that, at those energies, the most ef- 
ficient thermalization of neutrons takes place during the 
first interactions with the polyethylene moderator. Con- 
versely, the three rings converge to nearly the same effi- 
ciency at energies above 1 MeV, where the total efficiency 
drops significantly (see Fig. [7]). 

In the present experiment, the energies of the /3- 
delayed neutrons E^n ranged from zero to Qp — Sn, where 
is the /3-decay Q value of the mother nucleus and Sn 
is the neutron separation energy of the daughter nucleus. 
The distribution of Epn between these two values follows 
^ Sf3{E)f{Qp - E), where SpiE) is the /3-decay strength 
function for a decay into the daughter's level at energy 
E = Sn + Ei3n, and 



f{Qp ~E)^ (Qp -Sn- Epnf 



(5) 



The strong energy dependence of / largely favors Sp to 
excited levels of the daughter nucleus near its Sn- More- 
over, as discussed in Refs. [13, [U, H^l, high-resolution 
spectroscopic studies of /3-delayed neutron-emitter nu- 
clei produced by fission showed that Efjn was always 
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much lower than — Sn (e.g. 199 keV for ®^Br, 
450 keV for ^^Rh, and 579 keV for ^^h). This trend 
was also observed by the same authors in the total spec- 
tra of ^^^U and ^^^Pu, with average Ejjn of 575 keV 
and 525 keV, respectively, and with little neutron in- 
tensities at E(3n>S00 keV [sEIIi. The reason for these 
"compressed" Ep^ spectra is the strong, often preferred 
population of the lowest excited states in the final nu- 
clei [82| . Since the region investigated in the present 
work includes strongly deformed nuclei, the respective 
expected low-laying excited states are rather low. Thus, 
it was safe to assume the E/Sn values of the nuclei of inter- 
est to be typically below 500 keV. For these low energies, 
a constant value of (37±5)% for the NERO eSiciency was 
assumed. 



2. Neutron background 

Free neutron background rates were independently 
recorded throughout the experiment using NERO in self- 
triggering mode. Four of these measurements were taken 
without beam, and two with beam on target. The back- 
ground rates doubled from 4 to 8 s~^, when fragments 
were sent into the experimental setup, revealing the ex- 
istence of two different background sources. The energy 
spectra recorded during the background measurements 
proved that the background origin could be attributed to 
actual neutrons. One of the neutron-background sources 
was intrinsic to the detector and its environment, while 
the other had a beam-linked origin. Analysis of the ring- 
counting ratios for background and production runs sup- 
ported this idea (see Fig. [5]). Measurements of /9-delayed 
neutrons emitted from the implanted nuclei showed that 
the NERO counting rates were higher for the innermost 
ring (i.e., the closest to the DSSD) and systematically de- 
creased for the next external rings. This results is com- 
patible with MCNP simulations summarized in Fig. [T] 
Background runs with beam off showed the opposite 
trend, with high rates in the most external ring, which 
gradually decreased for the next internal ones. Such a 
result suggests that these runs were mainly affected by 
an external background source, most probably related 
to cosmic rays. Background runs with beam on target 
showed an intermediate situation that could be explained 
as arising from a combination of external and internal 
sources. 

The value of in Eq. [1] included contributions from 
neutron- /3-background events (i.e., neutrons in coinci- 
dence with /3-decay background events) B^^Bp), and 
from random coincidences between free NERO back- 
ground events and real (3 decays B^{f3). The value of 
BniBp) for each nucleus was calculated as the product 
of the neutron-/3-background rate measured on each im- 
planted DSSD-cluster, and the neutron-detection time 
following the corresponding implantations of that nu- 
cleus. Owing to the very low total number of neu- 
tron and /3-decay background coincidences measured per 



DSSD-cluster, the neutron-/5-background rate on each 
cluster was determined by scaling the /3-decay back- 
ground rate calculated in Sec. lIII A"T1 The corresponding 
scaling factor, calculated as the DSSD cluster-averaged 
ratio of neutron-/3-background coincidences to /3-decay 
background events, was about 0.08 and nearly constant 
throughout the experiment. Besides this background 
source, i?n(/3) was approximately calculated as the prod- 
uct of the number of mother /3-decays Np, and the prob- 
ability for at least one free neutron background with a 
rate 8 s~^ to be detected in random coincidence with 
a [3 decay. This latter approximation is not valid for 
coincidences of /3-decays with free background neutrons 
that were produced by fragmentation reactions induced 
by the same implanted mother nuclei. A calculated prob- 
ability for this scenario, however, demonstrated that the 
occurrence of such a type of coincidences was negligible. 
Table ini shows the value of iV/j, A^„, B^{Bp) and B„(/3). 



3. Error analysis 

The error analysis of P„ was derived from Eq. 3) In 
general, the main source came from uncertainties in the 
number of detected /3-delayed neutrons Np-o^ and back- 
ground events, with typical values about 20% for each. 
The former had statistical origin, whereas the latter 
was calculated from the /3-background uncertainties, de- 
scribed in Sec. lIII A 11 and the error in the determination 
of the 0.08 scaling factor described in the previous sec- 
tion. An additional contribution of 15% to the total error 
came from uncertainties in the number of mother decays 
Np, which were calculated by propagating the uncertain- 
ties in e^, according to Sec. IIII A 31 Finally, an average 
13.5% relative error in e„ was calculated as described in 

Ref. [H,!!!]. 

In the case of ^"^Mo and ^^°Mo — where contributions 
from the daughter nuclei to the total number of /3-delayed 
neutrons was significant — the systematic error was gov- 
erned by uncertainties in the value of C in Eq. 21 The 
latter was derived from the error propagation of all the 
variables in Eq. [31 being Pun the main contribution. Rel- 
ative uncertainties of about 46% and 35% were obtained 
for ^°^Mo and ^^"Mo, respectively. 

The Pn values and their errors obtained in the present 
experiment are listed in table [TTl and systematically pre- 
sented for each isotopic chain in Fig. [51 Pn values were 
only deduced for nuclei with a statistically significant 
number of detected neutrons, i.e., with a number of de- 
tected neutrons above the number of background neu- 
trons plus uncertainties within a 1-cr confidence level. 
Otherwise, only upper limits of Pn were deduced using 
the method described in Ref. [s^l — for a Poisson distribu- 
tion of detected neutrons — with the extension proposed 
by Hebbeker to include systematic uncertainties in the 
input quantities [8^1 (see vertical lines in Fig. [5]). The 
upper limits were calculated for a confidence level of 68%. 
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FIG. 8: Ratio of neutrons detected with different NERO rings for three different runs: production (left), background with 
beam off (center), and background with beam on (right). Histogram bin numbers 1, 2 and 3 correspond to ring ratios R2/R1, 
R3/R2, and R3/R1, where 7?i_3 are the innermost, intermediate and external rings (see text for details). 
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FIG. 9: Measured Pn-values (full circles) and established upper limits (segments) for Y, Zr, Nb, Mo and Tc isotopes, compared 
with results from previous experiments [3^. [3^. IgiI] (empty circles). For the sake of clarity, the latter were shifted to the right 
by 0.15 units. The data are compared with two versions of the QRPA model of MoUer et ai: the version described in Ref. [411 ] 
(solid line), and the interim version QRPA06 described in Sec. UV] (dashed line). (See text for details). 



IV. RESULTS AND DISCUSSION 



The newly measured /3-decay half-lives included in 
Tableland Fig. [6] (as full circles) follow a systematic de- 
creasing trend with neutron-richness, agreeing, for most 
of the cases, with previous measured values (empty cir- 
cles). In some particular cases, half- lives from /3-decay 



isomers were found in the literature. In particular, Khan 
et al. [59] reported two different half-lives for ^""^Yei, 
presumably from low- and high-spin /3-decaying isomers. 
The Ti/2 measured in the present experiment for this 
nucleus is compared in Fig. [6] with the value found by 
Wohn et al. [60] , presently assumed to correspond to the 
ground state |6l|. Similarly, in the case of ^"^Yga, two 
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TABLE II: Pn values obtained in the present experiment. The resuhs are compared with available data from previous experi- 
ments (Literature), and with the versions QRPA03 and QRPA06 of MoUer's QRPA model. (See text for details). 



Isotope 


Nf3 




B^{Bp) 


Bn(/3) 




Pn 


(%) 














Present exp. 


Literature 


y^txrPMo 41 




lOOy 


58 


1 


0.4 


0.1 


< 10 


0.92(8) [61] 


0.2 


0.3 


lOly 


231 


3 


0.6 


0.4 


< 4 


1.94(18) [61J 


0.6 


1.3 


loay 


373 


10 


1.1 


0.6 


6(2) 


4.9(12) [61J 


3.4 


1.6 


lOSy 


185 


6 


0.5 


0.3 


8(2) 


8(3) [32] 


4.0 


4.1 


104y 


40 


5 


0.1 


0.1 


34(10) 




4.8 


3.9 


105y 


8 


1 


0.03 


0.01 


< 82 




22.4 


17.0 


1032r 


856 


5 


2.4 


1.4 


< 1 




0.0 


0.0 


1042, 


1470 


10 


4.3 


2.4 


< 1 




0.0 


0.0 


i052r 


529 


4 


1.5 


0.8 


< 2 




0.0 


0.0 


106 2r 


199 


4 


0.5 


0.3 


< 7 




0.7 


0.7 


lo^Zr 


28 


1 


0.04 


0.04 


< 23 




0.6 


2.1 


106Nb 


3238 


70 


9.1 


5.2 


5(1) 


4.5(3) [32] 


0.3 


0.2 


lo^Nb 


2068 


68 


5.7 


3.3 


8(1) 


6(2) [32] 


4.4 


3.7 




458 


15 


1.3 


0.7 


8(2) 


6.2(5) [32] 


15.6 


11.0 


109Nb 


83 


3 


0.2 


0.1 


< 15 


31(5) [32] 


13.6 


26.0 


i°«Mo 


5557 


35 


24.2 


8.9 


< 0.5 




0.0 


0.0 


i«9Mo 


2856 


27 


8.1 


4.6 


1.3(6) 




0.0 


0.0 


"»Mo 


689 


8 


1.9 


1.1 


2.0(7) 




0.0 


0.0 


"iMo 


52 


1 


0.1 


0.1 


< 12 




0.0 


0.1 


109tc 


906 


6 


2.6 


1.4 


< 1 


0.08(2) [32] 


0.0 


0.0 


"OTc 


2960 


14 


8.5 


4.7 


< 4 


0.04(2) [32] 


0.2 


0.1 


"iTc 


1684 


12 


4.7 


2.7 


< 1 


0.85(20) [32] 


0.4 


0.6 




371 


6 


0.5 


0.6 


4(1) 


1.5(5) [33] 


0.9 


0.8 



half-lives were separately reported for the low-spin [63 
and high-spin [gs'] isomers. Interestingly enough, only the 
latter case is compatible with the value measured in the 
present experiment, thus indicating a favored production 
of this nucleus in a high-spin configuration. 

Our results include new half-lives for the iV=66 mid- 
shell isotopes "'^"^Ygg and "'^"^Zree, as well as the more 
exotic ^'^^Zrgy and -^^-^Moeg. New Pn values were also de- 
duced for "^°*Y65 and ^°^'-^^°Mo67,68j and new upper hm- 
its for lo^Yee, ^^^'^^^ liQz-m and i°«'"iMo66,69. In the 
case of ^^'^Ygs , the evolution of Pn with neutron number 
shows a pronounced increase compared with the smooth 
trend observed for lighter isotopes. Conversely, the sharp 
increase of Pn observed by Mehren and collaborators |33 l 
from Pn=(6.2±0.5)%, for lo^Nber, to Pn=(31±5)% for, 
^''^Nbgg, is not supported by our measured upper limit 
Pn<15% for lo^Nbgg. 

The small Pn values for ^"^'^^"Mogj^esi which could 
not be observed in previous experiments, were detected 



here as a result of a lower neutron background rate 
of 0.001 . Additionally, the selectivity achieved in 
the present experiment, resulting from the combined 
in-flight separation technique and the event-by-event 
implantation-decay-neutron correlations, made it possi- 
ble to rule out any potential neutron contaminant from 
neutron emitters in the cocktail beam. On the contrary, 
Wang et al. [sl] pointed out the possible presence of neu- 
tron emitting contaminants from neighboring isobars in 
IGISOL-type experiments, which can only be detected 
by continuous monitoring of 7-lines from the separated 
beam. These authors use that argument as a possible ex- 
planation for the weak components of long-lived contam- 
inants in the v4=104 time-spectrum measured by Mehren 
et al. [11]. 
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QRPA results 

The experimental data shown in Fig. [SJ for and 
in Fig.[9l for Pn, are compared to two theoretical calcula- 
tions. The solid lines are the results taken from Ref. [41} 
(QRPA03), in which the allowed Gamow- Teller transition 
rates are calculated in a microscopic QRPA approach and 
in which the first-forbidden transition rates are obtained 
from the statistical gross theory [s^, [s^- The second 
calculation (henceforth referred to as QRPA06), repre- 
sented by the dashed lines, gives results from an identi- 
cal model, but with the theoretical Qp values which enter 
in the phase-space integrals obtained from the irnproved 
finite-range liquid-drop model (FRLDM) of Ref. |43|| cor- 
responding to the last line in Table I of Ref. [43]. In 
this interim global mass model, triaxial deformation of 
the nuclear ground state is taken into account, but there 
are also some other improvements. The agreement with 
the 2003 mass evaluation is 0.6038 MeV. As noted else- 
where [3] , there are substantial effects from axial asym- 
metry on the ground-state masses in precisely the re- 
gion of nuclei studied in this work. When experimental 
masses are available for both parent and daughter Qp 
and Sn are calculated from experimental data, otherwise 
from theory. In QRPA03, the 1995 mass evaluation from 
Ref. jm was used, in the QRPA06-interim calculation, 
the 2003 evaluation [s^ was used. Finally, both mod- 
els calculate Sp assuming the same deformation for the 
mother and daughter nuclei; an approximation that was 
discussed in detail by Krumlinde and MoUer for some se- 
lected cases 38 1 . For complete details about the model 
si, HO, Ei • Examples of how different types 



see Refs. m, 



of nuclear structure effects manifest themselves in the 
calculated Ti and Pn are discussed in detail in, for ex- 
ample, Refs. [33, [M m [n m, m, m, m ) . 



General trends 

As shown in Figs. [5] and [51 QRPA06 shows gener- 
ally better agreement with the measured T1/2 and Pn 
than the older QRPA03. The generally poor results for 
the half-lives of the less exotic isotopes are consistent 
with the fact that uncertainties in parameters such as 
Qp have a very strong impact for decays with small en- 
ergy releases. This general behavior was already observed 
and discussed by Pfeiffer et al. for different nuclei (see 
Figs. 8 and 9, and Tables V and VI of Ref. [27|). Be- 
yond these general observations, the level of agreement 
between measured data and calculations shows no clear 
general systematic behavior. For instance, the half-lives 
predicted by both models are too short for all Y isotopes, 
and too long for all Mo isotopes. A similar trend is seen 
within the same isotopic chain such as Mo, where the 
half-life of ^"^Moee is well reproduced, while the more 
exotic "'^°^Mo67, ^^"Moes and ^^^Mogg are significantly 
overestimated. Such "fluctuating" behavior stems from 
the wide variety of nuclear shapes in this shape-transition 



region. 

Both QRPA03 and QRPA06 predict half-lives for all 
A^=65 isotones that are too short relative to the observed 
data, as was already pointed out by Wang et al. in their 
analysis of ^^'^Ygs [33|- According to these authors, the 
coupling of the proton orbital 7r[422]5/2+ to the neutron 
valence orbital j/[413]5/2"'" — which is in near proximity to 
j^[532]5/2^ at quadrupole deformation £2— 0.3 — would 
give rise to the allowed /3-decay transition from ^"'^Ygs 
0^ into the ^°^Zr64 0"*" ground state with a very short 
half-life. This interpretation explains also the disagree- 
ment between our measured and calculated Pn for iV=65 
isotones. In this case, the too low Pn values predicted 
by QRPA reflect an overestimated /9-decay feeding into 
levels below the neutron separation energy S'n. 



Analysts of nuclear deformations 

In the case of (3 decay of the Y isotopes, both mod- 
els behave similarly, showing too short half-lives and too 
low Pn values. The improved treatment of deformation 
in QRPA06 had no major impact when compared with 
QRPA03, as triaxiality is not expected to develop for 
these nuclei. Indeed, spectroscopic studies of Zr isotopes 
between A''=60 and A^=64 showed that these nuclei are 
dominated by increasing prolate deformations with no 
indication of triaxial components [sO, [9ll, [H, In 
an attempt to extend the analysis of nuclear shapes be- 
yond ^°^Zr64, we have re-calculated the T1/2 and Pn 
of ^°^Y65 and ^^^Yge, assuming different pure prolate 
shapes for the corresponding mother-daughter systems. 
Results from this analysis are shown in Fig.[Tni where the 
measured ri/2 and Pn are compared with calculations 
performed over a large range of quadrupole deformation 
(—0.35 < e2 < 0.35) of the daughter nuclei. 

Three remarks from this analysis, in regards to the (3 
decay of ^'^''Ygs into ^°^Zr64. Firstly, the calculated val- 
ues of Ti /2 and Pn experience an abrupt transition from 
their maxima, for a spherical daughter ^'^^Zr64, to very 
low values at deformations around e2— 0.1. Secondly, ex- 
perimental ri/2 and Pn are reproduced assuming a pro- 
late deformation e2'^0.20. Thirdly, for larger deforma- 
tions beyond e2>0.25 the (3 decay becomes faster with 
decreasing probabilities for /3-delayed neutron emissions. 
The good agreement of the calculations at e2~0.20 is 
ruled by GT transitions into four-quasiparticle levels at 
energies around S'n. Conversely, the too low predicted 
Ti/2 and Pn are governed by the fragmentation of Sfj 
into low-energy (i.e., well below S'n) two-quasiparticle 
states involving the coupling of 7rgg/2 levels with t^Si/2 
(at €2—0.1) or with high-fi Nilsson orbitals (at 

£2^0.25). Finally, the high Pn values around 100% for /3 
decay into a spherical ^°'*Zre4 arises from one single GT 
transition to the TrgQf2 ® ^^57/2 level at 7.25 MeV, well 
above Sn. 

Similarly to ^^-^Yes, the /3-decay half-life of ^"^Yee 
can only be correctly reproduced for deformations of the 
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FIG. 10: QRPA-calculated /3-decay half-lives and Pn values as a function of the quadrupole deformation of the daughter nuclei 
£2 (solid line), compared with measured data within 1-a uncertainty (shaded area). 



daughter ^"^^Zres given by e2~0.2, whereas the Pn upper 
limit serves only to rule out spherical deformations. At 
£2—0.2, the calculated T1/2 is governed by transitions into 
prolate three-quasiparticle 7r(7g/2(8> 2^^5/2 states at energies 
around S^- Conversely, at larger and smaller deforma- 
tions, the calculated /3 decay is too fast due to the pres- 
ence of different low-lying one-quasiparticle states that 
push Sf3 down in energy, well below Sn- Although the 
measured T1/2 is also compatible with oblate deforma- 
tions £2^—0.3, we ruled out such scenario as there is no 
experimental evidence of oblate shapes in lighter Zr iso- 
topes [9S|9ll,[93,[9|. 

In summary, the measured T1/2 and Pn of ^"''Yes 
and ^^^Yge can only be correctly reproduced if one as- 
sumes €2^0.2 for the corresponding daughter ^'^*Zr64 
and ^°^Zr65. On the basis of the larg er elongations 
found for lighter Zr isotopes [o^, [sHj l92j | . the new val- 
ues of £2 for ^°'*Zr and ^^^Zres contradict somehow the 
rather expected maximum saturated deformations at the 
iV=66 mid-shell. Interestingly enough, the maximum al- 
lowed quadrupole deformation of ^^^Ztq4 deduced from 
our data (e2~0.25) disagrees with the large P2 —0.4 



value obtained from analysis of the quadrupole moment 
Qq of the yrast band [9^, and from measurements of 
Since Sfs is sensitive to the nuclear 
structure of the daughter nucleus — including, besides the 
yrast band, any other level from its ground state to en- 
ergies just below Qi3 — our result points to the possible 
presence of spherical or weakly-deformed low-lying bands 
coexisting with a highly-deformed yrast band of ^''^Zr64. 
Spectroscopic studies of ^^"Zreo by Mach et al. re- 
vealed such coexisting bands with /32'^0.4 and /32'^0.2. In 
addition, analysis of the Qo and B{E2) systematics for Zr 
isotopes from iV=50 to A^=64 by Urban et al. [9l[ showed 
that coexisting spherical or weakly-deformed structures 
may be present beyond ^"^^Zrgo, although these authors 
claim that such phenomenon may end at iV=64. Since 
the QRPA formalism used in our analysis does not in- 
clude deformation of excited levels, the e2~0.2 presented 
here should be considered as an "effective" ground-state 
deformation resulting from the mixture of weakly and 
highly deformed bands in the daughter nuclei. Such a re- 
sult suggests that shape coexistence may still be present 
at i°'^Zr64 and ^"^Zrgs. This in turn may reflect the "tail- 
ing effect" of the predicted re-occurrence of the Z=40 
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sub-shell, together with a new sub-shell N=70 very far 
from stability [13, [H, , or the development of a more 
exotic tetrahedral shape at ^^"Zryo [29^ . 

The calculated T1/2 and for ^^^Zrge, ^"^ZrgT and 
^^^Mo69 as a function of 62 for the corresponding mother- 
daughter systems are shown in Fig.[Tni Here again, com- 
parisons between the measured T1/2 with calculations as- 
suming pure quadrupole deformation allow to constrain 
the possible €2 values of ^"^Nbgs and ^"''Nbge, whereas 
the calculated Pn values show not enough variation to 
distinguish between different deformations within uncer- 
tainties. The calculated T1/2 for ^^^Moes disagrees with 
the data for any pure quadrupole deformation of ^^^Tcgs- 
In this context, the FRLDM model predicts a triax- 
ial component 7=15° for ^o^Nbes, ^"^Nbge and "^Tcgs, 
which agrees with the values deduced by Luo et al. for 
^°^mei (7=13° [13) and "^Tcgs (7=26° (H). Interest- 
ingly, the measured T1/2 and Pn of ^^^Zrge and ^'^''Zrer 
are in excellent agreement with the results from QRPA06 
(which includes triaxiality) , as shown in Fig.[Bl Similarly, 
the calculated ri/2 of "'^"'^^Moeg is significantly improved 
when triaxiality is included, although no agreement with 
the measured value was yet found. 

In summary, analysis of the new measured data for 
lo^Zrge, ^°^Zr67 and "^Mogg using our QRPA06 calcula- 
tions indicates triaxial deformations for the correspond- 
ing daughter nuclei ^"^Nbes, ^^^Nbge and "^Tcgs. 



V. SUMMARY 

We have reported on the measurements of /9-decay 
properties of neutron-rich Y, Zr, Nb, Mo and Tc, which 
include new half-lives for ^°^Y, ^oeAor^j-^ and "^Mo, 
along with new P„ values for lo^Y and i^^'^^^Mo and P„ 
upper hmits for i03-i072j. and i°8^"iMo. The new data 
could be attained due to the low /3-decay background and 
/3-delayed neutron background rates obtained with the 
BCS/NERO detection setup. The high selectivity of the 
A1900 in-flight separator at NSCL was also an instru- 
mental achievement for the unambiguous identification 
of the exotic nuclei, thus demonstrating the optimum ca- 
pabilities of this experiment setup to reach very exotic 
regions, near (and at) the r-process path. 

The half-lives were analyzed using the MLH method 
and, in cases with enough statistics, least-squares fits 
of the decay curves. Agreement between both analy- 
sis brings confidence to the results. Analysis of the 
measured ri/2 and Pn based on QRPA model calcula- 
tions brings new insights to explore this interesting re- 
gion in terms of deformations. The measured ri/2 and 
Pn of ^°^'^''^Y65,66 isotopcs could Only be reproduced for 
quadrupole deformation parameters €2 of the correspond- 
ing daughter nuclei ^*''*'^°'^Zr64,65 below the values re- 
ported in the literature for ^°'*Zr64 and lighter isotopes. 
Since the /3-strength function Sjs governing the /?-decay 
from ^°^Y65 and ^"^Ygg is sensitive to the level struc- 
ture of the corresponding daughter nuclei, we believe 



that the low 62 derived in the present work for Zrg4 
and ^°^Zr65 is a probable signature of coexisting weakly- 
deformed bands. Such an interpretation is supported 
by previous independent analysis of deformations based 
on measurements of yrast-band quadrupole moments Qq 
and B{E2) for Zr isotopes between A^=50 and N—6A. 
The deformations reported in the present paper, however, 
show that weakly-deformed bands may still be present 
for -^°'*'^°^Zr64,65- The persistence of shape-coexistence 
for ^^'^ZrQi and ^^'^Zrgs may indicate the existence of a 
(near-) spherical doubly-magic ^^"^Zryo nucleus; a result 
that is compatible with the quenching of the iV=82 shell 
gap necessary to correct the unreahstic A~110 r-process 
abundance trough predicted by r-process model calcula- 
tions. 

The QRPA calculations also show that triaxial shapes 
play a critical role in the /3-decay to ^"^'^"^Nbgs^ee and 
^^^Tcgs. The inclusion of this new deformation degree 
of freedom — on the basis of the new FRLDM of MoUer 
et al. — significantly improves the calculated ri/2 and Pn 
with respect to the new measured values. In addition, 
the FRLDM-predicted triaxial components are compat- 
ible with values reported in the literature for nuclei in 
this region. 

Extension of /3-decay and spectroscopic experimental 
studies to full r-process nuclei requires new high-intensity 
fragmentation-beam facilities like FRIB and FAIR. These 
measurements are necessary to understand the nuclear 
physics governing the r-process. New measurements of 
masses and Qfj values of N—82 r-process isotones below 
Sn will clarify the role of shell quenching in the synthesis 
of heavy nuclei. 
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APPENDIX A: MAXIMUM LIKELIHOOD 
ANALYSIS OF HALF-LIVES 

The maximum likelihood method is the mathematical 
correct description even in cases of poor statistics. Let 
a set of independently measured quantities Xi originate 
from a probability density function f(xi,a), where a is 
a set of unknown parameters. The maximum likelihood 
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method consists of finding a parameter set a which max- 
imizes the joint probabiUty density 

i 

for aU measured data points Xi. C is also called the 
likelihood function. In most cases, it is easier to use InC 
instead and to solve the likelihood equation 



of simplicity, we will exclude in the present discussion the 
/3-delayed neutron branchings. 

Let Ai, A2, and A3 be the decay constants for the de- 
cay of a mother, daughter and granddaughter nucleus, 
respectively. It is important to distinguish the probabil- 
ity for a decay within a time t characterized by a decay 
constant Ai 

Fi(Ai,i) = l-e-^i*, 



d\nC 
da 



0. 



Normalization factors, which depend on the set of pa- 
rameters OL have to be included in the maximization pro- 
cess. All other multiplicative constants in the f{xi,a) 
can be neglected, even if they depend on the measured 
quantities Xi. 

The individual decay events of a decay sequence are not 
statistically independent, therefore the likelihood func- 
tion has to be defined. Additional corrections must be 
used to compensate for the neglected late decay events, 
if the correlation time window is small compared to the 
mean life time of the investigated nuclei. A method only 
using the first measured decay event within the correla- 
tion window is reported in [73] . Such a method does not 
make use of all available information and therefore might 
be disadvantageous in the case of poor statistics. Based 
on the work in [99] , the mathematical correct probability 
density function for up to three decay events within the 
correlation time window was developed. 

We assumed that all the decay events after an implan- 
tation of an identified nucleus within a position and time 
correlation interval belong to the first three decay genera- 
tion (mother, daughter and granddaughter decay). Addi- 
tionally, coincidentally assigned background events with 
a constant rate might occur. Up to three events within 
the time correlation window are considered. For the sake 



from the probability density function for a decay at 
exact time t characterized by a decay constant Ai 



/i(Ai,i) = Aie 



-Alt 



For the detection of the second decay generation, we 
use the probability for a decay within a time t of a daugh- 
ter nuclei with a decay constant A2 , which was populated 
by a mother decay with decay constant Ai: 



F2(Ai,A2,i) = 1 



A1A2 
A2 — Ai 



-Aii 



Ai 



A2 



and the probability density function for a decay of a 
daughter nuclei with a decay constant A2 at time t, which 
was populated by a mother decay with decay constant Ai : 



/2(Al,A2,i) = 



A1A2 



A2 — Ai 



3-Ait 



-,-A2t 



Similarly, the probability F3, and the corresponding 
probability density function /s, for a decay within a time 
t of a granddaughter nuclei with a decay constant A3, 
which was populated by a mother and daughter decay 
characterized by decay constants Ai and A2 are given by: 



F3{Xl,X2,X3,t) = 1 



A1A2A3 



(A2-Ai)(A3-Ai)(A3-A2) 



(A3-A2) ^-Ait _ (A3-Ai) ^_A2i ^ 

Ai A2 



(A2 — Ai) Aa* 



As 



and 



/3(Al, A2, A3,t) 



Ai A2A3 



(A2-Ai)(A3-Ai)(A3-A2) 



[(A3-A2)e-^i* - (A3-Ai)e-^^* + 



+(A2-Ai)e" 



Finally, for background events, the average rate and time is known. The probability for the observation of 
the expected number of events within the correlation exact r background events within a correlation time tc 
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and a background rate b can be calculated using Poisson 
statistics: 

= • 

Depending on the number of observed decay events 
within the correlation time, one has to consider aU possi- 
ble scenarios leading to the observation. In the following, 
we use a short notation to identify the composition of the 



probability terms of the various scenarios. Di stands for 
the probability that a decay of the i-th generation occurs, 
Oi that an occurring decay is observed, ei , €2 and €3 des- 
ignate the detection efficiencies for the respective decays, 
this is the probability for the observation of an occurring 
decay. In addition, the notations F{X, t) = 1 — F{X, t) 
and e = 1 — e are used. 

The probability for the observation of no decay event 
within the correlation time can be calculated as follows: 



Po(Ai) = (Di + + D1O1D2O2DZ + D1O1D2O2DMB0 

Po(Ai) = [A(Ai,tc) + (^'2(Ai,A2,ic) -^'i(Ai,ic))ei + (^'3(Ai,A2,A3,<c)- 
— F2(Ai, A2, tc))eie2 + P3(Ai, A2, A3,tc)eie2e3]Bo 
= [1 - Fi(Ai,ic)ei - F2(Ai, A2,ic)eie2 - i^3(Ai, A2, A3,ic)eie2e3]Bo. 



For the case of the observation of only one decay event 
within the correlation time, four scenarios are possible: 

1) The decay of the mother was observed, daughter 
and granddaughter decay did either not occur or these 
decays were not observed: 

PlOl = P{dl) = DiOl{D2 + D2O2D3 + D202D303)Bq. 

2) The decay of mother and daughter did occur, but 
only the daughter decay was observed, whereas the 
granddaughter decay did not occur or was not observed: 

P102 = P{d2) = Di6iD202{D3 + D303)Bo. 



3) All three decays did occur, but only the granddaugh- 
ter decay was observed: 

P103 = P{d3) = D1O1D2O2D3O3B0. 

Until now, we assumed that there was no background 
event within the correlation time. 

4) The last scenario describes the observation of a 
background event, all three decays did not occur or were 
not observed: 



Pi04 = P{b) = {Di + D1O1D2 + D1O1D2O2D3 + DiOiD202D303)Bi. 



r 



Calculation of the likelihood function requires the single decay event at time ti: 
probability density functions for the observations of a 
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Pioi(Ai) = Ci/i(Ai, ti)ei [Fi(A2, - h) + (F2(A2, A3, - h) ~ Fi{\2,t, - h)) 

£2 + -F2(A2, A3, tc — ti)e2l'i\ Bo 

= Ci/i(Ai, ti)ei [1 - Fi(A2, t, - ti)e2 - F2(A2, A3, - h)e2e3] Bo, 

Pi02(Ai) = Ci/2(Ai, A2,ii)eie2[A(A3,tc - h) + Fi(A3,tc - h)e3]Bo 
= Ci/2(Ai, A2, ii)eie2 [l - i^i(A3, - h)e3]Bo, 

Pi03(Ai) = Ci/3(Ai, A2, A3, ii)eie2e3-Bo, 

Pl04(Al) = Ci [A(Ai, + (F2(Al, A2, ie) - Fi(Ai, + (i^3(Ai, A2, A3, ie)- 

- i^2(Ai, A2, t,))eie2 + F3(Ai, A2, A3, t,)eie2e3] Sii^^ 
= Ci [1 - Fi(Ai, tc)ei - F2(Ai, A2, ic)eie2 - -FalAi, A2, A3, ic)eie2e3] 



The joint probability density function for observing one 
decay event at time ti is the sum of the single probability 
densities: 

Pi(Ai) = pioi(Ai) +pio2(Ai) +pio3(Ai) +pio4(Ai), 
where the normalization constant C'l fulfills the equation: 

/ pi(Ai)dti = l. 
Jo 

I 



Ten different scenarios need to be considered when 
two decay events occur within tc, and twenty different 
scenarios for three decay events. A detailed descrip- 
tion of all scenarios and the resulting normalized joint 
probability functions p2(Ai) and p3(Ai) can be found in 
Refs. [99. 100]. 

The analysis program assigns decay events to preceding 
implantations, only events with one (n^ = 1), two (rii — 
2) or three = 3) decay events within the correlation 
time are considered. Therefore, we initially maximize the 
likelihood function for A''i23 observed decay sequences: 



JV123 

/:i23(Ai) - n ['^("^ ~ l)Pi(Ai) + 6in, - 2)p2(Ai) + <5(n, - 3)p3(Ai) . 
1=1 

The solution Ai^, of the maximization equation 



dXi 



= 



has to be corrected for events with no observed decay 
events within the correlation time window. The most 
likely number of events A^o with no observed decay events 
within the correlation time depends on Po(Ai) and there- 
fore on Ai itself: 



Nn = 



PoiXi) 
l-Po(Ai 



23- 



To find the solution Xij_^-^^ of the maximization equation 
of the joint likelihood function C, an iterative numerical 
method is used until Ai converges: 



/:(,+i)(Ai) = £i23(Ai)Po(Ai) 



i9/3(j+i)(Ai 



dXi 



0. 



Ai — Ai 



The correlation time, therefore, should be long com- 
pared to the mean life time of the mother nuclei to avoid 
large correction factors due to this iterative method. If 
the background rate is low enough, a correlation time 
equal to ten half-lives should be used. If the time win- 
dow is too long, the assumption of a maximum number 
of three decays within the correlation time is no longer 
valid and the maximum likelihood method might fail. 
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The validity of this approximation as well as a check of 
the whole procedure for the analysis of decay sequences 
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